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We study numerically the topological knot solution suggested recently in the Weinberg-Salam 
model. Applying the SU (2) gauge invariant Abelian projection we demonstrate that the restricted 
part of the Weinberg-Salam Lagrangian containing the interaction of the neutral boson with the 
Higgs scalar can be reduced to the Ginzburg-Landau model with the hidden SU (2) symmetry. The 
energy of the knot composed from the neutral boson and Higgs field has been evaluated by using 
the variational method with a modified Ward ansatz. The obtained numerical value is 39 Tev which 
provides the upper bound on the electroweak knot energy. 
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1. Introduction 

It has been well known since 1980's that knot-like 
topological solitons corresponding to a non-trivial Hopf 
number exist in some non-linear er-models with a scalar 
field n a ,a — 1,2,3 UI 12, |3|- Last decade a wide class of 
such topological knot solitons have been proposed in var- 
ious physical modcls: in the effective theory of quantum 
chromodynamics (QCD) m two-gap superconduc- 

tivity models [jj, and in two-component Bose-Einstein 
condensate theory 0, H, @- Recently it was suggested 
that knot solitons could exist even in high energy physics, 
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namely, in the Weinberg-Salam model 

In the present paper we consider numerically the 
properties of the knot soliton made of a neutral Z bo- 
son with a dressing Higgs field in the Weinberg-Salam 
model. We apply the Ward ansatz which provided a 
high accuracy for the energy of the Faddeev-Niemi knot 
and results in, as we will show below, a reasonable value 
39 Tev for the energy of the electroweak knot. Even 
though the obtained energy value seems to be out of 
practical detection with the present experimental facil- 
ities, the electroweak knot may have an impaet on cos- 
mological problems related to dark matter and baryon 
asymmetry of the universe 

[H IH IH- In particular, it 
was suggested that the formation of primordial hyper- 
magnetic knots and their further decay could result in 
the baryon asymmetry 0, 0. We suppose that the 
existcncc of the topological knot in the Standard model 
can lead to an alternative mechanism of the origin of the 
baryon asymmetry in the universe, like the baryogenesis 
seenario based on electroweak strings 0] . 

2. Gauge invariant Abelian projection 

First we consider the SU (2) Yang- Mills theory with 
a complex Higgs scalar doublet as a simple model of 
Weinberg-Salam theory. We will extract from the La- 
grangian of the Weinberg-Salam model the part contain- 



ing the interaction of the neutral boson with the Higgs 
scalar using the gauge invariant Abelian projection. By 
this way we kccp cxplicitly the full SU (2) gauge invari- 
ance and at the same time we preserve the topological 
structure of the original group SU (2) in the Abelianized 
part of the gauge connection. This approach provides the 
gauge invariant definition of the Hopf number in terms 
of the Abelian potential and gives rise to the deseription 
of the neutral boson interaeting with the Higgs scalar in 
the framework of a simple Ginzburg-Landau model. 

We will follow the approach proposed in [Tol Il7| . 
Let us start from the gauge invariant Abelian projection 
which is given by the following decomposition of the full 
SU (2) connection [l7| 
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The restricted potential A^ possesses transformation 
properties of the full SU (2) connection and leads to the 
desired Abelianization of the theory. The veetor part of 

the restricted potential C M = — h x d^h is defined in 

terms of the triplet scalar field n whose configuration is 
classified by Hopf invariant, i. e. by a homotopy group 
7T3(5 2 ). The scalar field n is covariantly constant 

D^h = d^n + gA^ Xn = (2) 

and represents the topological degrees of freedom in the 
theory. The SU (2) gauge transformations can be written 
as follows 
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Sn = —a x n , SAfj, — —D^a, 



(3) 



SAn = -h ■ d^a, 

g 



SAft = -D^a, 



-d x w u 



(4) 
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In a local orthonormal frame (ni, 77.2,^3 = n) of the 
internal SU(2) space one can rewrite the decomposition 
(jTJ as 

A, = A,n+{J^+W^)h {1 . 2) , (5) 

where w/j , : are local components of the vectors 
CpjW/j, correspondingly. By appropriate gauge transfor- 
mation one can pass to the magnetic gauge [l7| in which 
the vector potential C M converts into the Abelian mag- 
netic potential 



C> = --ni9 M n 2 . 
9 
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With the decomposition QJ one can factorize the 
Abelian potential part and the ofT-diagonal components 
in an explicit gauge invariant manner 

F„u = (Fpv + H MV )fi 

+£ M W>„ - D U W^ + gW^ x W v , (7) 



where 



F^u — d^A v d v A^, 
1 
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e abc n a d^n b d v n c = 3^-8^. (8) 



So that the Lagrangian can be rewritten in the form 
C = --F 2 

= -~(i> + H^f - \{b»w v - b v w»f - 

2 

+ H MV ) ■ (W» x W„) - ^-(W^ x W,) 2 . (9) 

In complex notation for the charged vector boson 
1 



(10) 



one can express the Lagrangian explicitly in terms of the 
Abelian potential A^ + C M and the charged vector field 

C = —{F^ + H^f - \\D^W V - D V W^\ 2 

2 



-^DpWv-DvWn] 2 , 



(H) 



where 



Wp, = -ig(Wn W u— W v 
DyWu = {d^+igiA^ + C^Wu. 



It should be stressed, that the Lagrangian is invariant 
not only under the original SU(2) group transformations 
but under the action of the dual U(l) group as well 



s su(2)(A f _ l + C M ) 



0. 



(12) 



Due to this the Hopf number can be defined as a Chern- 
Simons index of SU (2) gauge theory. For the case A^ = 
cmployed later on we have the Standard formula for the 
Hopf number |l7| 
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Notice, the Abelian symmetry U(l) has been appeared 
in QCD as a dual magnetic symmetry [l^ which is an 
essential ingredient of the dual Meissner effect leading to 
the monopole mechanism of confmcmcnt. 

With these preliminaries let us start from the bosonic 
sector of the Weinberg-Salam Lagrangian 
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where is the Abelian potential corresponding to the 
weak hypercharge group E/V(l). 

We consider the restricted case with the vanishing 
charged vector boson — and electromagnetic field 
A^ = while keeping only the Higgs field <p and the 
neutral Z boson in the Lagrangian. We will identify the 
classical part of Z boson containing the whole topological 
content of the knot with a hypermagnetic potential C M 
|l0j . It is convenient to express the fields h, <fi in terms 
of the real components p, a and SU (2) complex doublet 
£ which parametrizes the complex projective space CP 1 



e + e = 1. 



(15) 



Since the hypermagnetic potential C„ is the dual poten- 
tial one can express it locally as follows 
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(16) 



With this we obtain the final expression for the La- 
grangian for our purpose to study the static knot type 
solution 
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(17) 



3 



Formally, the expression l|17(l can be easily recog- 
nized as a well-known Ginsburg-Landau Lagrangian. But 
one should stress that the original non-Abelian struc- 
ture of the Lagrangian (|17f) is hidden and encoded in 
the topological properties of the hypermagnetic poten- 
tial C p, given by Ijl6|) . By this the non-trivial homotopy 
group tt^CP 1 ) ensures the existence of the topological 
knot solution classified by Hopf number N u ■ 

3. Knot soliton 

The idea of a non-topological st ring in the electroweak 
theory was proposed long time ago . Later the twisted 
electroweak strings with a baryon number had been in- 
troduced in [l^, and the existence of the knot in the 
Weinberg-Salam model was suggested recently [lC|| . We 
will find the static topological knot solution by mini- 
mization of the energy functional corresponding to the 
Lagrangian l(T7|) (i = 1,2,3) 



for the complex doublet £ (in spherical coordinates) 
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where h(r) = f (r) /r. With this one can write down the 
components of the magnetic potential 



C f — 
C e = 



2 h' sin 2 9 cos 9 
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2h sin 3 9 
r(l + h 2 ){h 2 +cos 2 e)' 
2 sin 9 
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(23) 
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Let us introduce the dimensionless variables r — > 
rmw, k = 2%/A/g = mn/mw, Ci — > Cig/mw, p — > 
p/v with rescaling a — > «cos^vk & n d tow — 80.22 Gev , 
mff = 200 Gev, sin 2 8w — 0.2325. One can rewrite the 
energy functional in terms of the dimensionless variables 
as follows 
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where = tij^djC^ is the hypermagnetic field 

strength of the neutral boson. 

To find out appropriate trial funetions we start from 
the original Ward ansatz The stereographic projec- 
tion of n is given by a complex field 
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ni + in 2 



l + n 3 

The Ward ansatz with a radial trial funetion f (r) 



(20) 



z - i f (r) ' 
f(r)=a(r-b)(r 2 
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was proposed to deseribe the Faddeev-Niemi knot soli- 
ton with a good accuracy for the energy value. The non- 
trivial Hopf number N u — 1 is determined by the bound- 
ary conditions of the funetion f (r) at points r = 0, oo. 
Applying this ansatz with the trial funetion f (r) and the 
definitions (|15f) one can find the corresponding expression 



One can use the gauge freedom to remove the depen- 
dence of the phase faetor ip(r, 8, p) on the variables (r, 9). 
By making an appropriate gauge transformation 



Ci -> C t - 2Vi arctan( -), 

cost^ 

ip(r, 9, p) -> ip(r, 9, p) - arctan(— ^— ) (24) 

cos 9 

one results in a simple parametrization for the phase fae- 
tor ip(r, 9, p) — p and for the components C r , Cg 



2ti cos 9 
1 + h 2 ' 
2hsin9 

r(l + h 2 )' 



(25) 



The phase faetor in the scalar field <j> — pe la appears 
in the Lagrangian only in the combination Ci + 2<9i<x 
We will adopt the following ansatz for the phase fae- 
tor a(r,6) = (3 (r) cos 9 with the boundary conditions 
/3(0) = /3(oo) = 0. The chosen form of the ansatz pre- 
serves the Hopf number caleulated with the gauge trans- 
formed potentials 
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It turns out that the modified ansatz l|26|l provides lower 
value of the minimum of the energy functional. For the 
Higgs field p(r, 9) we impose the boundary conditions 



p(0,9) = c, p(oo, 



= 1, 



(27) 



where c is a number. The minimization of the energy 
functional gives the energy of the knot Eq = 39 Tev. This 



FIG. 1: The contour plot of the energy density of the knot in 
the plane y — 0. The maximal value of the energy density is 
£max = 7.42 , the step size between the contours is 0.39 (in 
units of v/g = 0.407 Tev/m 3 ). 



energy value has been checked with a class of various trial 
functions including the angle dependence on 9. We found 
that the minimum of the energy functional is reached for 
a sphcrically symmetric configuration of the Higgs field. 
We expect that our modified Ward ansatz H2bj) with a 
pure radial trial funetion for the Higgs field p(r) provides 
the main contribution to the knot energy similar to the 
case of the Faddeev-Niemi knot. The energy density plots 
in the plane XY and XZ are shown in Figs. 1, 2 from 
which one can see that the energy of the knot is localized 
along the torus. 

The topological stability of the knot is provided by 
non-vanishing hypermagnetic fluxes along and around 
the torus. The total hypermagnetic fiux around the axis 
0Z is ^ = 2ir. The total hypermagnetic flux through 
the inner disk of the torus (r < tq = 1.64,0 = tt/2) is 
<3> 2 = 2tt. By this we have the linking of two fiuxes which 
provides the Hopf number Qh = 1 and the topological 
stability. Notice, that the lines of the hypermagnetic vec- 
tor field H form in faet a twisted torus. 

One should stress, since we use a variational method 
with a constrained ansatz we can not rather obtain a 
more detailed structure of the energy density distribu- 
tion of the knot. By plotting the energy density one 
can deseribe only qualitatively the inner structure of 
the knot and estimate roughly the effective size of the 
knot. The energy density has a maximum along the cir- 
cle ro = 0.50, 8 = tt/2 which corresponds to the radius of 
the knot 1.45 x 10 _18 m. From the other hand we can de- 
fine the radial size of the knot as a radius of the cirele at 
which the hypermagnetic component Hg changes its sign. 
This gives the radius r = 4.1 x 10~ 18 m which is near 
the value 3.5 x 10~ 18 m estimated in [Io|. Notice, that 



FIG. 2: The energy density £(x, y, 0) of the knot in the plane 

z = 




FIG. 3: (a) the energy density £{r,ix/2) (in unit of v /g); 
(b) the hypermagnetic field H g (r, 7i"/2); (c) the hypermagnetic 
field H v (r, n/2), (d) the Higgs field p(r); (e) the phase radial 
funetion /3(r) (in dimensionless units). 



the Ward ansatz gives exactly the Dirac type magnetic 
flux quantization with the minimal non-zero flux <f> = 2ir 
due to the semi-simple structure of the electroweak gauge 
group SU (2) x and the non-compaetness of the 

electric charge group U e (l). For a model with a sim- 
ple gauge group SU (2) we would have Wu-Yang type 
flux quantization with a flux quantum proportional to 
47r. Such situation has been realized, for instance, in 
a two-component Bose-Einstein condensate system |9j. 
The energy density and hypermagnetic field components 
Hg,H v in the plane 6 = tt/2, the Higgs field p (r) and the 
phase radial funetion (3 (r) are shown in Fig. 3. 

Since the energy of the knot is much greater than 
the energy of veetor W, Z bosons an important question 
arises whether the knot solution is stable against the in- 
teraetion with other fields, especially with the charged 
W boson and electromagnetic field. The only poten- 
tially dangerous part in the Weinberg-Salam Lagrangian 
which can cause the instability is represented by the 
so called anomalous magnetic moment interaetion term 

igH^v W 'n W v in (fTTf) - This term leads to a severe vac- 
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uum instab ility problem in QCD via generating a tachy- 
onic mode [2(j . 

To analyse the stability of the knot soliton we consider 
small fluctuations of W field in knot background. The 
possible unstable modes are determined by the negative 
eigenvalues of the operator K defined by means of the 
second variational derivative of the Lagrangian with re- 
spect to infinitesimal variations SW^. In a suitable back- 
ground gauge (d^+ig cos ^CJ^ = one can find the 
operator K as follows (in dimensionless variables) 



K lut = - 



,9 Z 



V^{ d p d p - P 2 ) + i-p z tw 



(28) 



where, D,, 



tCu. One can factorize the eigen- 



values of the matrix part with Lorentzian indices in the 
operator K similarly to the case of derivation of func- 
tional determinants corresponding to the one-loop effec- 
tive action in S U (2) QCD (details of derivation for gen- 
eral magnetic and electric background are given in |2l|) 



K* ~ -D t D t + p'± 2-H, 

g 



(29) 



where H = \H\. With C M defined by l|26[) one can write 
down the eigenfunetion equation for the operator 



\d r (r 2 d r ) - —L—d 9 {smOde) - 

J.Z r i g m gl 

c 2 



9' 



did 



9 



12 



c 



, 2e 
p 2 ±-H 

g 



^(r,0) = A* (r, 6») (30) 



The equation represents a Schrodinger type problem for a 
quantum mechanical partiele moving in a potential well. 
To estimate the lowest eigenvalue corresponding to the 
ground state one can neglect the terms 2iCidi and idiC'i 
which do not contribute to the ground state in the lead- 
ing order approximation due to the axial symmetry and 
linear dependeney of the knot solution C, on angle func- 
tions. The negative eigenvalues can arise only from the 
equation with the lower negative sign in front of the last 
term in l.h.s of (|30[l . With this, substituting the mi- 
merical knot solution into the equation and solving it by 
the variational method one can find the lowest eigenvalue 



which is found to be a positive number Ao = +0.31 . The 
scalar field p plays a role of the mass regulator parameter 
which is elose to the mass of W boson (see Fig.3). One 
can find numerically that negative modes would appear if 
p were replaced with a constant averaged value less than 
0.73. Surprisingly, despite on the faets that the knot en- 
ergy is much higher than the energy of W boson, and our 
solution is obtained with a simple modified Ward ansatz, 
nevertheless, no any negative eigenvalues appear. The 
weak dependeney on the knot energy is provided by the 

e 2 

presence of two competetive terms in (I30II . —p^Cf and 



2e 



g'- 



H, which contribute opposite ways to the quantum 

9 . 

mechanical potential in the Schrodinger equation. Wc 
have checked also that unstable modes do not appear 
under small deviation of the profile funetions /, p, a from 
the numerical solution. This analysis gives us a hope 
that a rigorous exact solution for the knot will be stable 
as well. 

In conclusion, we performed numerical study of the 
topological knot soliton in the Weinberg-Salam model 
suggested in [ljj. The obtained numerical value of the 
knot energy is 39 Tev which is higher than the earlier 
estimate 21 Tev ^(J. The last estimate was based on 
the equation for the knot energy obtained in the non- 
linear a-model Q which represents actually the low en- 
ergy bound. Since we elaborate the energy of the knot 
with the variational method, our numerical result pro- 
vides the upper bound which, we believe, is elose to the 
real value for the energy of the electroweak knot. The 
knot energy value even being much higher than the scale 
100 Gev of the Standard model, it is still less than the 
scale of the temperature T ~ 80 Tev for the right-electron 
cquilibration state considered in estimation of the baryon 
asymmetry [2^, [23j. This might be an indication that 
the topological electroweak knot could play an impor- 
tant role in cosmology at the stage of formation of large 
scale structures along with other proposed mechanisms. 
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